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INTRODUCTION 


The motion of a sphere free falling in the high atmosphere deviates from gravita- 
tional acceleration according to the drag expression 

PV^C^A 

where a^ and g^ are the inertial and gravitational acceleration components along 
the trajectory. (Symbols are defined in the appendix.) Atmospheric density can be 
deduced from the motion according to 


_ 2m ~ 6y 
CdA 

Uncertainty in the result depends on uncertainty in the measured trajectory quantities 
assumed drag coefficient, gravitattonai acceieration, and sphere constants. By assum’ing 
random independence, these uncertainties are related as follows: 



Atmospheric temperature can be calculated from the deduced density profile within 
the validity of the hydrostatic equation and equation of state according to 


T = 


Mw 


P 



Uncertainty in the initial value Tq and mean molecular weight, as well as in the 
aforementioned quantities, produce uncertainty in the calculated temperature. 
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Where according to the nonindependence of the error In g with altitude 

and correlation in such quantities. 

ACCELERATION AND VELOCITY UNCERTAINTY 

uncertainty in the measured acceleration and velocity of a passive sphere is deter 
■ a h Le accuracv of the tracking instrument and by the length of the data-smoothing 
mrrval used in computing time derivatives. (We make the tacit 
is introduced by the smoothing model, which is tantamount to assuming ttet th p 

Hvr nf thp order or form of the polynomial or other function used in 

dItirBases for a general quantitative evaluation of 
from position uncertainty and smoothing interval lei^h are presen e 

The tracking system measures the position of the sphere to within some ™lume of 
1 V whL in turn may vary in slae and shape according to sphere position. The 

r s^^ronl^est Tn tMrl^ 

tX under discussion is parallel to the path tangent. By assuming that the tracking 
::::: measures sphere positi. in s^t .ng^ B aainiu„ 
angle E, with uncertainty in each, 6R, OA, ana or., u y 
path tangent is given by 

\2 


(5X)2 = 


( R5R)^ + (r^ESe) + (r^ASA cos e) 
V2 


where 


y2 _ + (RE) + (RA cos E) 


i.2 


Smoothing of data rests on available discriminating knowledge of the signal and 
■se L I data Thl maior characteristic of the signal, which immediately suggests 
Zu Ts sLeLe inertia, which Implies a cutoff frequency In sphere velocity response. 

Linearizing the drag acceleration about a point on the trajectory provides an expres- 

Sion for the local frequency response. 

pCdV 


Olg = 


(m/A) 


The corresponding minimum altitude wavelei^h of density structure to which 
can respond is 


Zfi = 


2irz 


Cilk 
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The radius 
given by 


Az of the altitude neighborhood within which a constant 


1 ( 

' dU)Q 

a)g' 

\ az 


^Az 


= + J_ 

\P 9z Cd 


8C 


D 


az 




ojg is valid is 


The first term in the sum is the reciprocal of the local scale height 


H = 

the second is small, and the third is proportional to density according to 

1^ W _ P^D V 
V az 2(m/A) z 

The first term dominates at high altitude. By assigning a "validity limit" Awg/wg, the 

altitude interval over which constant Wg holds is ±Az = The character of 

low-order polynomial smoothing routines in the frequency domain is approximated by that 
of a low-pass filter with cutoff frequency o)f inversely proportional to smoothing 
interval W. 

‘^f- W 

The altitude interval traversed during a "window" time length W is proportional to 
vertical velocity z, and the minimum altitude wavelength in density structure passed by 
such a filter is about twice that interval, 

Zf = 2Wz = ^ 


Processing the data with a low-pass filter having a cutoff frequency = u)^ 
would reject considerable noise without disturbing the signal. As altitude increases, 
however, Zg/2 exceeds its range of validity 2H(Act)g/o)g^ . Though approaches are con- 
ceivable, by using advanced processing techniques which may effectively enable the use 
of greater smoothing intervals without exceeding the validity of the associated smoothing 
model, smoothing intervals are currently limited for practical purposes to those over 
which, ujg is constant, i.e., to Zg/2 or 2H(Aa)g/wg), whichever is smaller. 


Matching W to the frequency response of the 
2H(Aojg/Wg^ requires 


Ty ^ 7r(m/A) 
~ pCdV 


sphere but limiting it to 

[Zg ^ 4H(Au)g/a;g)j 



z 


[^s > 4H(Awg/Wg^ 
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For the higher altitudes of interest, the window will depend on H, and therefore, the 
data are filtered less and less with increasing altitude, relative to sphere response. A 
variable smoothing interval of this maximum length approximates best that can be done 
with the standard smoothing techniques of low-order polynomial fitting. 

Other noise and signal characteristics which permit the identification of certain 
desired or undesired components in the data sometimes exist. Special techniques 
designed for specific components are assumed operative in the present discussion to the 
extent that their respective net noise components do not exceed the assumed measure- 
ment uncertainties. 

Time derivatives from a low-pass filter are bounded by the cutoff frequency and 
by input amplitude. If input uncertainty is 6X, then uncertainties in the first and second 
derivatives are approximated by 

6X = Wj6X 
6X = wjsx 


The above relations combine and reduce to an error model for the passive-sphere 

Ao>s 

system at high altitudes. By letting = 1, 



The first term is 


20 


(m/A) 


- 6X 


(Zf = Zs) 



f + 4 

\HV ) 

\2HpCd vJ 


(Zf = 4 h) 


A quantitative illustration of this error model, including the associated temperature 
uncertainty is presented in tables I and n and figure 1. 

For reference, a simulated 160-km trajectory in the U.S. Standard Atmosphere 
1962 of a 1 -meter, lOO-gram sphere is used. Arbitrary values for input uncertainties 
are based on published experience (refs. 1,2, and 3). 

Uncertainty is taken here to correspond to practical limits of error. Uncertainty 
of the value of a normally distributed random variable is taken as two standard deviations 
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and that of a uniformly distributed quantity its semirange (radius). Quoted uncertainty 
values for drag coefficient, for example, are understood to mean the limit within which 
the true value lies with a high (80- to 90-percent) probability. Inputs and results here 
may be halved throughout to correspond to the standard deviation. 

The uncertainty inputs are as follows; 

6R = 9,0 meters 


6E = 6A = 0.0001 radian 
5Cj^A 


CdA 


= 0.03 


gy 


= 0.01 


^= 0.01 

m 


6M 


M 


^= 0.03 


W 


OBSERVATION 

The quantitative results of the illustration suggest a significant finding concerning 
the passive-sphere technique. Deviations between sphere results and results from 
other atmospheric sounding techniques may be explained in large part as a consequence 
of the smoothing interval used. 

It is recommended that data reduction programs include the automatic calculation 
of the point-to-point uncertainty and present a measured parameter q from a given 
sounding in the form of q ± Aq, rather than simply as q. Thus, for example, the density 
profile plot would be a curve with varying width within which the measured density lies 
with a stated probability. Comparisons between techniques would then be expected to 
produce overlapping curves. 

It is expected that tracking accuracy, 6R, 6E, 6A, will not increase significantly 
over the assumed values for altitudes considerably higher than 90 km. Therefore, larger 
data-smoothing intervals would materially enhance the altitude capability of the passive - 
sphere technique. By assuming successful higher order smoothing models, an upper 
bound or the maximum possible smoothing interval at a given altitude is the period of the 
entire trajectory above that altitude including both upleg and downleg portions. Further 
investigation of the practical capabilities of potential and advanced data processing 
methods is indicated. The tracking accuracy requirements for a new passive -sphere 
system must depend on properties of the smoothing interval to be used. 
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APPENDIX 


A 

ay 

Cd 

E 

g 

gy 

H 

M 

Mw 

m 

NRe 

R 

R* 

T 

t 

V 

W 


SYMBOLS 

sphere cross-sectional area, m^; radar azimuth angle, radians 
drag acceleration, m/sec^ 

inertial acceleration along the trajectory, m/sec^ 
drag coefficients 

radar elevation angle , radians or degrees 
gravity, m/sec^ 

gravitational acceleration along the trajectory, m/sec2 

atmospheric scale height, m 

Mach number 

mean molecular weight 

sphere mass, 

Reynolds number 
radar slant range, m 

universal gas constant, joules/*^K kg-mol 

atmospheric temperature, °K 

time 

air speed, m/sec 
smoothing interval, sec 
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X position along trajectory, m 

Zf,Zg minimum altitude wavelength; passed by smoothing or sensed by 

sphere, respectively, m 

z altitude, m 

6q uncertainty along trajectory of any parameter q 

fig degree of nonindependence of gravity error with altitude 

p atmospheric density, l^/m^ 

\)/ trajectory angle from horizontal, degrees 

ojjjOjg cutoff frequency of smoothing response or of sphere response, respectively, 

radians/sec 

Subscripts: 
i,j,k indices 

o initial 

Dots over a symbol indicate the degree of the derivative with respect to time. 
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TABLE I 


DOWNLEG 


t, 

sec 

km 

mjs’ec 

'i', 

degrees 

R, 

km 

degrees 

m/sec 

RE, 

m/sec 

m/sec 

220 

115 

960 

-72.8 

143.9 

53.0 

-562 

-779 

-917 

225 

110 

1004 

-73.3 

140.9 

51.3 

-575 

-821 

-963 

230 

105 

1043 

-74.2 

138.0 

49.6 

-583 

-856 

-1004 

235 

100 

1074 

-74.8 

135.1 

47.8 

-582 

-911 

-1037 

240 

95 

1089 

-75.9 

132.4 

45.4 

-569 

-925 

-1055 

245 

90 

1073 

-75.7 

129.7 

43.9 

-539 

-930 

-1042 

250 

85 

986 

-76.9 

127.2 

42.0 

-475 

-862 

-960 

255 

80 

774 

-77.4 

124.8 

39.5 

-358 

-718 

-756 

264 

75 

458 

-79.2 

122.4 

37.8 

-208 

-406 

-450 

279 

70 

254 

-83.0 

120.1 

35.7 

-123 

-221 

-252 


km 

M 

Nr6 


m/?ec2 

115 

2.75 

2.41 

2.44 

0.413 

110 

3.12 

5.99 

2.37 

.914 

105 

3.40 

14.5 

2.20 

1.98 

100 

3.69 

38.4 

1.84 

4.13 

95 

3.88 

101.0 

1.51 

8.52 

90 

3.98 

277.0 

1.29 

18.4 

85 

3.66 

645.0 

1.19 

36.1 

80 

2.87 

1277.0 

1.12 

53.1 

75 

1,61 

1491.0 

1,03 

36.8 

70 

.854 

1540.0 

.604 

13.3 


2, 

km 

p 

6T 

T 

6T, 

OK 

6X, 

m 

6aj3 

2^ 

V 

Zf, 

km 

2sj 

km 

w, 

sec 

115 


0,909 

275,0 

12.8 

0.757 

0.004 

36.9 

6. +6 

20.1 

no 


.593 

153.0 

12.6 

.519 

.005 

31.2 

3. +6 

16.3 

105 

.314 

.349 

81.6 

12,5 

.313 

.005 

28.4 

1.+6 

14.1 

100 

.199 

.214 

45,0 

12.4 

.197 

.006 

25.4 

847 

12.2 

95 

.117 

.126 

24.0 

12.2 

.113 

.006 

23.6 

424 

12.8 

90 

.067 

.076 

13,7 

12.1 

.059 

.007 

21.8 

191 

10.5 

85 

.041 

.052 

9.4 

12.0 

.025 

.007 

21.8 

82.4 

11.5 

80 

.034 

.046 

8.3 

11.8 

.011 

.007 

21.7 

34.6 

14.3 

75 

.033 

.046 

9.2 

11.7 

.008 

,008 

24.0 

17.6 

19.2 

70 

.034 

.046 

10.1 

11.4 

.009 

.009 

26.3 

15.1 

30.0 
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TABLE II 


UPLEG 



km 

m/sec 

degrees 

R, 

km 

E, 

degrees 

m/sec 

RE, 

m/sec 

m/sec 


115 

962 

72.8 

119.0 

75.0 

961.7 

-35.9 

919 


110 

1012 

73.8 

113.9 

75.0 

1012.0 

-24.7 

972 

13.4 

105 

1065 

74.4 

108.7 

75.1 

1065.0 

-13.8 

1026 

8.6 

100 

1124 

74.9 

103.5 

75.1 

1124.0 

-2.32 

1085 

4.1 

95 

1197 

75.5 

98.2 

75.2 

1197.0 

8.93 

1159 

0 

90 

1304 

76.0 

93.1 

j 

75.2 

1303.0 

20.6 

1265 


z, 

km 

M 

NRe 

Cd 

m/sec^ 

115 

2.76 

2.40 

2.44 

0.414 

110 

3.15 

6.04 

2.37 

.930 

105 

3.47 

14.8 

2.19 

,413 

100 

3.87 

40.5 

1.82 

.442 

95 

4.27 

112.0 

1.48 

10.2 

90 

4.84 

340.0 

1.27 

26.8 


z, 

km 

P 

6T 

T 

6T, 

OK 

6X, 

m 

6ap 

2^ 

V 

Zf, 

km 

Zs> 

km 

W, 

sec 

115 

0.535 

0.735 

223.0 

9.00 

0.534 

0.003 

36.9 

6. +6 

mm 

110 

.371 

.468 

120.0 

9.00 

.369 

.003 

31.3 

3.+6 

16.1 

105 

.228 

.273 

63.9 

9.00 

.226 

.004 

28.4 

1.+6 

13.9 


.148 

.166 

34.9 

9.00 

.144 

.004 

25.4 

8.+5 

11.7 

95 

.090 

.100 

19.6 

9.00 

.084 

.005 

23.6 

427 


90 


.065 

11.7 

9.00 

.045 

.005 

21.8 

193 

8.6 
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